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I consider the non-equilibrium transport of electrons through a quantum system with a thermo- 
electric response. This system may be any nanostructure or molecule modeled by the nonlinear 
scattering theory which includes Hartree-like electrostatic interactions exactly and certain dynamic 
interaction effects (decoherence and relaxation) phenomenologically. I derive three fundamental 
bounds for such quantum systems coupled to multiple macroscopic reservoirs, one of which may be 
superconducting. These bounds affect nonlinear heating effects (such as Joule heating), and work 
and entropy production. Two bounds correspond to the first law and second law of thermodynamics 
in classical physics. The third bound is quantum (wavelength dependent), and is as important as the 
thermodynamic ones in limiting the capabilities of mesoscopic heat-engines and refrigerators. The 
quantum bound also leads to Nernst's unattainability principle that the quantum system cannot 
cool a reservoir to absolute zero in a finite time, although it can get exponentially close. 

PACS numbers: 73.63.-b, 05.70.Ln, 72.15.Jf, 84.60.Rb 



I. INTRODUCTION 

In classical physics, currents in resistive circuits are 
dissipative; they radiate Joule heat which increases the 
entropy of the environment of the circuit. Devices 
with thermoelectric responseiP^ can be used as heat- 
engines (using heat-flows to create electrical power) or 
heat pumps (using electrical power to create heat-flows) . 
These too are dissipative, since they increase the entropy 
of the environment (at least a little bit) during their oper- 
ation. The laws of thermodynamics give us strict bounds 
on the capabilities of such thermoelectric devices. 

There is currently great interest in using quantum sys- 
tems as thermoelectric he at-engines or refrigerator, par- 
ticularly nanostrutures^HHI and molecules^^. How- 




arbitrary nanostructure 




arbitrary molecule 

Figure 1: Quantum systems with K n — 3 normal reservoirs 
and one superconducting (sc) reservoir. Quantum systems 
may be (a) nanostructures such as those defined by electro- 
static gates on top of a two-dimensional electron gas (2deg), 
or (b) molecules. Charge and heat current, U and Ji, are 
positive if they flow from the quantum system to reservoir i. 



ever quantum responses are often non-local, since ther- 
mal equilibration occurs on lengthscales very much big- 
ger than the size of the system, with equilibrium only 
being established deep in the reservoirs coupled to the 
system. Classical models of transport account neither 
for this non-locality, nor for the wave-nature of quantum 
particles. It is natural to ask what the bounds are for 
these quantum devices, and how they compare to the 
bounds on classical devices. There are many works ask- 
ing whether quantum systems are bounded by the laws of 
thermodynamics^^, particularly for systems modeled 
by Markovian or Bloch-Redfield master equations, see 
Rcfs. [20 28-30J and references therein. However the fun- 
damental bounds on thermoelectric quantum systems re- 
mains an open question^. Here, I consider the fully non- 
linear response of a broad class of thermoelectric quan- 
tum systems, and show that they obey the bounds im- 
posed by the laws of thermodynamics, while also suffering 
a quantum bound on their behavior. 

The Landauer-Biittiker scattering theory^lMI j s a 
widely-used theory for modeling such electron flows 
through qu antum systems; giving DC and AC 
conduct ance d 35 * 3 ^ , and nonlinear effect:P3EED n has 
equally been used to get linear thermal and thermoelec- 
tric effects^H^. However it is a delicate question whether 
it captures dissipation and the associated laws of thermo- 
dynamics (or their quantum equivalent). The standard 
answer is that the theory accurately models reversible 
scattering processes, while assuming that equilibriation 
occurs (irreversibly) in the reservoirs via dissipative pro- 
cesses which the theory does not detail. Hence, one might 
suspect that the theory would not capture how dissi- 
pation produces entropy. This suspicion is wrong, the 
nonlinear scattering theory is sufficient to give informa- 
tion about entropy production, despite the simplistic way 
in which it treats the physics in the reservoirs. Indeed, 
Ref. [27] showed the second law of thermodynamics to be 
a remarkably direct outcome of a nonlinear "DC" scatter- 
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ing theory^ for any quantum system of non-interacting 
electrons between two reservoirs 32 . Energy conservation 
leds to the systems also obeying the first law. 

Here I show that the proof in Ref. [27 works equally 
well for a non-equilibrium scattering theory which in- 
cludes Hartree-like electron- electron interactions (electro- 
static charging effects) in a self-consistent and gauge- 
invariant manner. This theory is the analogue for heat- 
currents^ of the nonlinear scattering theory for charge- 
currents^Eni although heat-currents are not conserved 
when charge-currents are. Such electrostatic interactions 
cannot be ignored for the nonlinear response^^l, and 
entropy production is a purely nonlinear effect (vanishing 
in the linear-response regime). By keeping only electro- 
static interaction effects, this approach is by no means 
exact, however it is a reasonable approximation for any 
system where the single-electron charging energy is small 
compared to temperature or broadening of the quantum 
system's levels due to the coupling to the reservoirs^. 
Within this approach, I prove that the DC response of 
any quantum system coupled to an arbitrary number of 
reservoirs — one of which may be superconducting - 
obeys the laws of thermodynamics. 

The extension to many reservoirs (from two reservoirs 
in Ref. [27]) requires no new physics, but is technically 
complicated, because one can no longer use the simple 
form of two-terminal Onsager relations. However it ex- 
tends our physical understanding of thermodynamic laws 
in quantum-transport in two important ways. 

(a) Using the generalization to many reservoirs, one 
can include dynamics interaction effects (decoher- 
ence and relaxation) ph enomen ologically as addi- 
tional fictitious reservoir d 25 ! 47 * 48 !. Performing such 
an analysis, I conclude that the classical laws of 
thermodynamics hold for partially coherent (as well 
as fully coherent) quantum-transport. 

(b) A superconducting (sc) reservoir has no analogue 
in classical mechanics, because the condensate is a 
macroscopic quantum state. In addition, charge- 
current can flow into (or out of) this condensate, 
but heat-current cannot. Thus a charge-current can 
pass into the sc reservoir, / sc without any change in 
that reservoir's entropy; S sc — J sc /T sc = 0. Despite 
this, I show that the laws of thermodynamics hold 
in the presence of such a reservoir. 

I also derive a quantum bound (qb) on the thermo- 
electric response. It tells us that no matter how much 
DC power one supplied to any device, it cannot extract 
heat from the ith reservoir faster than the bound J 4 qb oc 
Ni x (ksTi) 2 , where Ti is the reservoir temperature, and 
the contact between the quantum device and the reservoir 
has Ni transverse modes. This bound is "quantum"; i.e. 
it is only relevant when A^ 1 ~ A|/Aj ^ 0, for electron 
wavelength Af, and contact cross-section A^. The quan- 
tity Jf h corresponds to the heat-current out of the zth 
reservoir when it is coupled to a zero-temperature reser- 
voir with the same chemical-potential, through a contact 



containing N modes. In such a context, it is a fermionic 
analogue^ of the Stefan-Boltzmann law of black-body 
radiation, for which J oc T . 

For thermoelectric refrigerators, the quantum and 
thermodynamic bounds compete, and the quantum 
bound can be a stronger constraint than that given by 
the thermodynamic bounds (the Carnot limit). In con- 
trast, the bounds combine to give an upper limit on the 
power that a thermoelectric heat-engine can produce. 

I find that the quantum bound means that the systems 
obey Nernst's unattainability principle (one of the two 
statements forming the third law of thermodynamics). 
This principle claims that a reservoir cannot be refrig- 
erated to absolute zero in a finite time, see for example 
Refs. |2ll29l30| for physical and Refs. |5D15T] for mathe- 
matical overviews. The quantum bound means that the 
systems considered here have (at best) the critical behav- 
ior; i.e. the reservoir cannot achieve zero temperature in a 
finite time, but it can get exponentially close. This rules 
out such systems as candidates for violating the unattain- 
ability principle in the manner Ref. [30_ proposes. 

II. EXPLICIT FORM OF THE BOUNDS 

In this section I give the three bounds, before deriv- 
ing them latter in the article. I consider an arbitrary 
quantum system coupled to any number of reservoirs, 
one of which may be superconducting. Normal reservoirs 
are treated as free-electron gases in thermal equilibrium, 
while a superconducting reservoir is treated as a con- 
densate of Cooper paira^l. The heat-current into reser- 
voir i is Ji. The electrical power flowing into reservoir 
i is Pi = Vili, with Vi being the bias on the reservoir, 
and Ii being the charge-current flowing into it. The rate 
of change of entropy in the ith reservoir is defined as 
Si = Ji/Ti for reservoir temperature Tj. 

I consider the DC response of the quantum system^, 
and assume that any variations of the energy or entropy 
of the quantum system over time are small enough to 
ignore. Defining, Stotal & n d <Stotai as the rate of total 
energy and entropy flow out of the quantum system into 
the reservoirs, averaged over long times. I show that 
scattering theory reproduces the thermodynamic laws, 

-Btotai - J2{Ji + Pi) = [1 st law] (1) 

i 

Stotal = V T < ^ [ 2 " d kW l ( 2 ) 

i 

where the sum is over all reservoirs. The quantum bound 
(qb) on the heat current into reservoir i is 

2 

Ji > -Jf = - Jr^Vi {Mf [quantum] (3) 
oh 

where the quantum system couples to reservoir i through 
N transverse modes, and this equation defines Jf h . Sec- 
tions |IX] and [X] discuss how these bounds act as limits on 
heat-engines and refrigerators. 
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III. NONLINEAR SCATTERING THEORY FOR 
HEAT TRANSPORT 



modulus squared of the elements of S^. . This property 
is that for all i, j, fi, v, e, v cnv , 



Here I consider the nonlinear version of the Landauer- 
Biittiker scattering theory of transport, which includes 
(in a gauge-invariant and self-consistent manner) those 
Hartree-type electron-electron interaction effects which 
can be considered as electrostatic charging effects. It has 
been used for AC charge-transporpSMl anc i nonlinear DC 
charge-transporlP^BQl Here, as in Ref. [13 J . I apply it to 
nonlinear DC heat-transport. 

One starttPmni 

with the scattering matrix, 
S^(e,v cnv ), of a ^-particle entering the quantum 
system from reservoir j to a /i-particle leaving via 
reservoir i, where \i and v can be either electron (with 
charge e") or hole (with charge — e"). The supercoduct- 
ing (sc) reservoir is modeled by Andreev reflection of 
electrons into holes and vice versa. Here an electron 
with energy e is an occupied conduction-band state at 
energy e + e~V sc , while a hole at energy e is an empty 
conduction-band state at energy — e + e~V sc , where 
e~V sc is the superconductor's chemical-potential. Thus 
these "holes" are not those in Ref. [T3], which are in a 
semiconductor's valence band. 

One assumes that the lead to reservoir i has iVj modes 
independent of e. The number of open modes in this 
lead can still be e-dependent; if a given mode n.j is closed 
off at a certain e, then it appears as purely reflected in 
the scattering matrix. If a superconducting reservoir is 
present, one measures the energy e from the supercon- 
ductor's chemical-potential e~V sc ] otherwise one can take 
e = however is convenient. I define v cnv as the set of 
all voltages and temperatures in the environment of the 
quantum system, {Vfc,Tfc}, where k runs over all reser- 
voirs and gates. In the nonlinear regime, S^(e, v env J 
must include the charging effects self-consistently, be- 
cause it is a function of the charge-distribution in the 
quantum system, which is in turn a function of § as we ll 
as the set of all voltages and temperatures in v en v ■ 

All the transport properties (both charge and heat 
currents), are then given by the transmission matrix, 

t oHV( 

y y-i ■ CI1V j j °ij 

the trace is over all transverse modes of the leads con- 
necting the quantum system to the i and j reservoirs. 
Any self-consistent analysis must of course respect gauge- 
invariance; i.e. shifting all voltages by the same amount 
is tantamount to redefining the zero of energy, and thus 
cannot change the physics. I ensure this by taking e 
and all reservoir and gate voltages (and any screening 
potentials which depend on them, as relative to the su- 
perconductor's chemical-potentialpSl e~V BC . If there is no 
superconductor, one can take these relative to another 
reservoir voltage. In everything that follows, I will as- 
sume that 7^ has this gauge-invariant. 

There are two properties of 7^ , which will be crucial 
in deriving the results presented in this article. The first 
property is due to the fact that 7^ is the sum of the 



7??( 



> 0. 



(4) 



7^(e v ) 



Tr 



[S« (e.Venv)] §f/(e,v cnv ) 



where 



The second property is due to the unitary of S^- (particle 
conservation), and is that for all e,v onv , 



A" 



j=l V 
= EE^Ve 



(5) 



i— 1 fj, 



where the i or j sums are over the K n normal reservoirs, 
while the /i or v sums are over ±1, with +1 for electrons 
(e) and —1 for holes (h). 

The charge current into the ith reservoir counts the 
charge (ie~ that enters or leaves that reservoir. The heat 
current into the ith reservoir counts the excess-energy 
(e — fie'Vi) that enters or leaves that reservoir-^. The 
charge current, Ii, and the heat-current, Ji, from the 
quantum system into the ith normal reservoir are 



■* = £ 



de 



h ^ 1 11 



M e-[T^(e)-^f/]/;( £ ), (6) 



de 



( e - Me -^)[^( e )-Ar^]//(e). (7) 



Here j is summed over the K n normal reservoirs, while 
fi and v are summed over e (+1) and h (—1). Here 8?? 
is shorthand for the product of Kronecker 8- functions, 
8ij8^ Ul while the Fermi function for j/-particles entering 
the quantum system from reservoir j, is 



//(e) = (1 + exp [(e - ueV 3 ) / {k B T 3 )\ ) 1 



(8) 



If a superconductor is present, then the charge-current 
into it (in the form of Cooper pairs), / sc , is given by 
Kirchoff 's law of current conservation. In contrast, heat 
current is not conserved. However heat flow through a 
boundary at the surface of a superconducting reservoir 
must be zero, since an electron at energy e is Andreev 
reflected as a hole with the same energy. Thus 



(9) 



For / sc the sum is over the K n normal reservoirs. 
Comparing Eqs. (|6| and Q, one easily see that 

J, = -V, h + J2j Q x e ^ (e) _ NiS ^ ^ (e) ' (10) 



Summing the heat current over all reservoirs (noting that 
J sc = 0), one can use Eq. (|5| to show that the sec- 
ond term in Eq. ( 10 1 cancels, so 53i=Ti Ji = ~ Y^i=i Vi 
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Note that I have chosen to measure all voltages relative 
to the superconductor (if present) . If instead, one mea- 
sured the voltages relative to an external reference, then 
the total heat-current out of the quantum system is 



'total 



i=l 



(11) 



Kirchoff's law in Eq. ^ means that changing all volt- 
ages by the same amount does not change the total heat 
flow, Jtotai- Thus Jtotai respects gauge-invariance when- 
ever J, does, where one recalls that Jj is a self-consistent 
function of 40 (Vi, T u V 2 , T 2 , ■ ■ ■ V N ,T N ). Whenever J total 
is non-zero, heat-currents are not conserved in the quan- 
tum system; the quantum system is a heat-source for 
Jtotai > 0, or a heat-sink for J tota i < 0. 



IV. ZEROTH LAW OF THERMODYNAMICS 

The definition of equilibrium is one of the various state- 
ments that together form the zeroth law. It is defined as 
a state for which one can cut the system into parts in 
any manner, and one will find no charge and heat flow 
between those parts. 

A consequence of this is that if a system in equilib- 
rium consists of multiple reservoirs (evidently in equilib- 
rium with each other) all coupled to each other through 
a quantum system, then there can not be any charge or 
heat current through that quantum system. It is trivial 
to use the scattering theory to show that any quantum 
system, placed between an arbitrary number of reservoirs 
all in equilibrium with each other, must obey this con- 
sequence of the zeroth law. It is sufficient to note that 
the scattering matrix § is unitary, this means that the 
transmission matrix obeys Eq. (|5| . Substituting this into 
Eqs. (6|7l, one see that charge and heat currents are zero 
whenever Tj — T Q and V 3 ■ — Vn for all j. 



V. 



FIRST LAW OF THERMODYNAMICS AND 
JOULE HEATING 



Eqs. ( 11 ) has the form of a classical Joule heating term; 
i.e. voltage x current. This is due to the fact that en- 
ergy conservation ensures that the heat-radiated (or ab- 
sorbed) is equal and opposite to electrical power supplied 
(or generated). Thus, one arrives directly at the first law 
of thermodynamics, Eq. 0. However it is important to 
note that this is only true upon summing over all reser- 
voirs, there is no particular relation between the power 
into reservoir i, Pi — Vili, and the heat-current into that 



reservoir, given in Eq. (10 I 



As a first example, consider a quantum system with 
two normal reservoirs at the same temperature but dif- 
ferent biases, Kirchoff's law means that I2 = —Ii, and 
hence Ptotai = (Vi — V<j)ii. Then p has the opposite sign 
to Vi — V2, and hence Ptotai < 0. This means the total 




rj"^)— |load2~[ — 



reservoirT"^)— | load 3 ~| — 
op)— [load sc I — 



sc reservoir 



Environment 
temp. = T a 



Figure 2: A quantum system with a finite thermoelectric 
response being used as a heat-engine, providing power to a 
set of loads, the sc reservoir is superconducting, but "load sc" 
is not (its simply the load coupled to the sc reservoir). 



heat current out of the quantum system, Jtotai = --Ptotai, 
is positive and corresponds to Joule heating. Note that 
if reservoir 1 is superconducting all the heat generated 
flows into reservoir 2. This means that superconductors 
can be used to filter out heat flow generated by Joule 
heating (as well as any other heat flow). 

As a second example, the quantum system in Fig. [2] is 
assumed to have a finite thermoelectric response. If one 
reservoir is heated it provides a total electrical power of 
Potai = Y2i Vj h to the loads. Then the total heat flow 
out of the quantum system equals —Ptotai (meaning more 
heat flows-in than flows-out) and the first law is satisfied. 



VI. SECOND LAW OF THERMODYNAMICS 

A. System coupled to two reservoirs 

Here I reproduce the calculation in Ref. [27] whose out- 
come is the second law in two-reservoir systems. The 
Hartree-type interaction effects can be included with- 
out any changes to their method. For two reservoirs, 
the rate of change of entropy is S = (J1/P1 + J2/J2), 
one can substitute in Eq. ^ and explicitly write out 
the sum over j — 1,2. There is no Andreev scattering 
so 7v, oc S^v, hence when there are only two-reservoirs 
Eq. ([5]) gives us the following two-terminal Onsager rela- 
tion, Tff = = 7$ = -7$. As a result, 



Stotal 



x T^(£,v cnv )[/(6 M )-/(^)L (12) 



where I define £ ifl = (e-^eVij/ik^Ti), so /f (e) = f(£ itl ) 
for /(£) = (1 + exp[£]) -1 . Since /(£) is a monotonically 
decaying function of £, the product of the two square- 
brackets in Eq. p2| cannot be negative. Taking this 
together with Eq. (Ul), one conclude that the second law, 
•Stotal ^ 0- is satisfied for such a two-reservoir system. 
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B. System coupled to any number of reservoirs 



1^1=1 l^i 



2K a v-^2if n 



m— 1 2-~d 



after which 



In general for more than two reservoirs, there is no 
such simple relationship between Tj/ 1 and as above. 
Then the derivation leading to the second law is more 
involved. Here I provide such a derivation; one that uses 



2K„ 



2K a 2K a 



only the unitarity of §f", in the form of Eqs. (45). Using 



Eq. Q, 



I write Stota! = (k B /h) J °° de Y(e) with 



Y(e) 



(13) 



where i and j sums are over the K n normal reservoirs. 
As in Eq. Q, I defined ^ = (e - fieVi) / (k B Ti) and 
/(£) = (1 + exp I now prove that Y(e) > for all 

e, thereby showing the second law, Eq. is satisfied. 
One can re-labeling the reservoirs in Eq. (13 1, giv- 



ing separate labels to the electron-state entering a given 
reservoir, and the hole states entering that same reser- 
voir. Thus I replace the sum E*=aEu=±i w ith E?=i > 
where i is related to i and /i as follows. For each i 
and p, one calculates and then order the elements 
in ascending order of One then labels these states 
with superscripts, i, running from 1 to 2K n , such that 
< £(2) < . . . < £(2**), i then define 

% 

£« = £(°)+^A£^, (14) 

n=l 

/( £») = /( £(0)) a/(™>, (15) 

m — 1 

where A£< s ) = [£< s ) - £(*-*)] and A/(™) = 
[/(£(™)) - /(£(™- 1 ))]. Here £<°) is arbitrary, so I take 
£(0) < £(i). ThenA£( K ) > for all n, and A/(" 1 ) < for 
all fh (since /(£) decays monotonically with £). 

Taking Eq. ( 13 I and replacing the quadruple sum over 
fi, v with a double sum oy er i, j, I substitute for £^ 
and using Eqs. ( |l4|l5| . Then 



2K„ 

He) = E 

*>5=l 



£W [ % (e)-JVi%] /(C (0) ) 
E A ^ /(C (0) ) 

n=l 

E C (0) [%(e) - Nify] A/(™) 

m— 1 

EE A ^ (S) [%(^)-^%] a/(™) 



n— 1 m— 1 



However given Eq. (Tsl) , the first three lines of this re- 
sult give zero when summed over i, j; leaving only the 
term containing the sum over both n and fh to evalu- 
ate. To do this, I write £~=i = J^lti E?5? and 



y( e ) 



£ A£(«) A/W E E 



n,m— 1 



z— n j—m 



For n > to, one can prove that the sum over j is negative 
(or zero) for all i. The proof requires noting that every 
sum over j contains a negative term for j — i and positive 
terms for j ^ i, and then using Eqs. ( 4|5 1 to show that the 
sum over J is not positive. This is not the case if n < fh, 
because not all sums over j contain a term with j = 1. 
However for n < fh one can do the i sum first; using the 
same logic as above to find that the sum over i is negative 
(or zero) for all j. Thus ^ (e) _ ^ < q 

for all n, to. Combining this with the fact that A£(") > 
for all n and A/(" 1 ) < for all to, tells us that Y(e) > 0, 
and thus the second law, Eq. ([2]), is obeyed. 

As an aside, I mention that the above proof that 
E?=5? Y.]=rn [%( e ) - N&j] < is most easily seen 
graphically. Write out the matrix [Tjj — NiSjj] explic- 
itly and draw a rectangle around the elements summed 
over for given n, fh. One sees that for n > m, every row 
of the rectangle contains a diagonal (z = j) element, so 
Eqs. (4j5) tell us that summing the elements in each row 
gives a negative (or zero) result. In contrast, ones sees 
that for n < fh, every column of the rectangle contains 
a diagonal (i = j) element, thus summing over the ele- 
ments in each column gives a negative (or zero) result. 



C. Positivity of Joule heating 

A very common situation is that where all reservoirs 
have the same temperature, so Tj = T for all i. Then 
the second law, Eq. reduces to the statement that 
the Joule heating radiated by any system, Eq. (Ill, will 
obey Jtotai > 0. 



VII. QUANTUM BOUND ON 
HEAT-CURRENTS 

Here I will take Eq. Q for the heat-current from the 
quantum system into reservoir i, and show that there is a 
negative lower-bound on its value. This shows that there 
is an upper-bound on the amount of heat that the DC 
thermoelectric response of a quantum system can suck 
from any given reservoir. This is irrespective of both the 
nature of the quantum system, and how much power one 
supplies to it. 

The proof relies upon splitting Eq. (|7| into two parts, 

(a) . . 

with J; L containing only those parts of the integrals 

where (e — fJ,e'Vi) > 0, and containing only those 
parts of the integrals where (e — [ieVi) < 0. In what 
follows I provide the proof for positive e~Vi, however the 
proof for negative e'Vi follows exactly the same logic (un- 
der the interchange e O h). For positive e'Vi, on e has 
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J L = J, (a) + j| b) where 



ji a) (e-^>0) 



• /e " v * j> 

/•oo J 



j| b) ( e -^>o) 



de 



where the ±1 is for e and h respectively. Our objective 
is to find the lower-bound on J. and . 

For J. (e"V< > 0), the first term in both integrands 
(the terms of the form (e±e~Vj)) is positive, while Eq. Q 
tells us that 7?" > 0, and the Fermi function satisfies 

< //(e) < 1- Thus jj b ^(e"Vi > 0) can never be more 
negative than if (e) = for all j ^ i and v ^ fi (where 
/i = +1 in the first integral and fx = — 1 in the second). 
The remaining terms — those with j = i and fj, = v — 
can never be more negative than if T^f = 0. Making the 
substitution £ = (e — fie'Vi) / '(ksTi) for /i = +1 in the 
first integral and —1 in the second, one arrives at 



4 a) (eVi>0) > 



N z (k B Ti) 



P(0) + P 



e-v; 



, (16) 



having defined P(x) = J^ d£ £/(l + exp[£]). 

For J^^e'Vi > 0), the first term in the integrand 
(the term of the form (e — e~Vi)) is negative. Thus 
J^ h \e~Vi > 0) is most negative when the j, z/-sum is as 
positive as possible, and this is when the Fermi func- 
tions fj(e) = 1 for all terms with j ^ i and v ^ +1. 
Then using Eq. ([5| to write the sum over (J,v) ^ (i, (i) 
in terms of the (z, /z) element of the sum, one sees that 

4 b) > I eVl f (e - erV0[^ - [l - /^W], 

which is most negative when ' + =0. Taking this 
lower bound, one can make the substitution £ = (e — 
e-Vi)/{k B Ti), and use the fact that [1 - /(£)] = /(-£), 
thus arriving at 



4 h \e-Vi>0) > 



P(0) - P 



eVj 



(17) 



Taking the sum of Eqs. ( |16|17 l, and noting that P(0) = 
— Lia(— 1) = 7r 2 /12, where Li2(z) is a dilogarithm func- 
tion, one arrives at Eq. Q3J) . Redoing the derivation for 

e"Vi < 0, gives different results for J> a ' and j| b ', but 
their sum is the same as for e'Vi > 0, thus the inequality 
in Eq. pi holds for all V z . 



VIII. PHENOMENOLOGICAL TREATMENT 
OF DECOHERENCE AND RELAXATION 

The scattering theory used above includes electrostatic 
interaction effects, but not dynamic interaction effects. 
This means that neither temperature nor bias induce the 
decoherence or relaxation (thermalization) of the elec- 
trons as they pass through the quantum system. How- 
ever, one can include these effects in the model in a phe- 
nomenological manner. The idea is to model a system 
with such dynamic interaction effects as follows. 

• Relaxation to a thermal state (and the associ- 
ated dec oherenc e) can be modeled by fictitious 
reservoir d 25 * 47 * 48 !, whose bias and temperature are 
chosen to ensure that the average charge and heat 
currents into each such reservoir is zero. 

• "Pure-dephasing" (decoherence without any energy 
relaxation) can be modeled by an infinite set of fic- 
titious reservoirs coupled to the system via leads 
which only lets a given energy pass through it, 
and tuning the bias on those reservoirs so that the 
charge and heat current into all of them are zero. 
The energy filtering in these leads can be treated 
as part of the system's scattering matrix. 

A system without dynamic interactions, but with the fic- 
titious reservoirs, obeys the bounds in Eqs. (l][3). Since 
the average charge and heat currents into the fictitious 
reservoirs are zero, the bounds on flows into the real 
reservoirs are the same as in the absence of the fictitious 
reservoirs. Thus it seems reasonable to conclude that 
Eqs. (T]|3l, are obeyed regardless of whether the electrons 
decohere and relax within the quantum system (due to 
electron-electron interactions) or not. 



IX. LIMITS ON QUANTUM HEAT-ENGINES 

Consider a two terminal thermoelectric device used to 
extract electrical power Ptotai from a heat flow between 
a hot reservoir (temperature T B ) and a cold reservoir 
(temperature 2c). The heat-current into the hot and 
cold reservoirs are Jh < and Jq > 0, while the power 
generated, Ptot ai > 0. Combining the first and second 
laws, Eqs. (T[2 1, one finds Carnot's result that 

Ptotai < - Jh (1 - Tc/Ph) ■ (18) 

This inequality limits the efficiency, rj = Ptotal/( — Jh); 
with which the heat-engine extracts electrical power from 
a given heat-current. The equality corresponds to Carnot 



efficiency. Combining Eq. (18) with the bound on Jh in 
Eq. (§, one gets the absolute limit on the power 

p„,., < ^W„- r c) (lg) 

For a heat-engine between a hot-source at Xh = 100°C= 
373K, and a cold source at T c = 20°C= 293K, the max- 
imum power per transverse mode is 15nW. 



7 



X. LIMITS ON QUANTUM REFRIGERATORS 

Consider a two terminal device extracting heat from a 
cold reservoir (to keep it cold or cool it further), when 
the ambient temperature is To and the cold reservoir's 
temperature is Tq < Tq. Power is supplied (negative 
-Ptotai) to extract heat from the cold reservoir (negative 
Jc). Combining the first and second laws, Eqs. (|1|2[ ), one 
gets the Carnot bound on refrigeration. Eq. (|3J) gives a 
second quantum bound. Together they read, 



--Ptotai T C /(T -Tc), [thermodyn. 



Jc < 



N c ir 2 (k B T c ) 2 /6h. 



[quantum] 



(20) 



The thermodynamic bound places no limit on the heat 
that one can extract from the cold reservoir; one can 
always extract heat if one is able to provide enough elec- 
trical power. In stark contrast the quantum bound places 
an absolute limit on the heat extracted. 

As a first example, consider a kitchen refrigerator run- 
ning at 100W, cooling from room-temperature Tq — 
293K to Tq = 0°C= 273K, and now imagine it consists 
of thermoelectric devices. Then the bounds are 



Jc < 



1300 W, 
A c x 36 nW. 



[thermodyn. 
[quantum] 



(21) 



The quantum bound is stronger than the thermodynamic 
one if Nq < 10 10 . A typical semiconductor thermoelec- 
tric (Fermi wavelength ~ lOOnm) with a cross-section of 
lcmxlcm, has Ac ~ 10 10 ; the quantum bound domi- 
nates for cross-sections smaller than this. 

As a second example, consider a point-contact nanos- 
tructure such as in Ref. |54| . There typical biases and 
currents are tens of mV and tenths of nA, so typical 
powers are pico- Watts (pW). Suppose such a device is 
used for cooling as proposed in Ref. [T3], to refriger- 
ate a micron-sized island from the cryostat temperature 
T - IK down to T c ~ 0.1K. Then 



J C < 



0.1 pW, 



[thermodyn. 



A c x 0.005 pW. [quantum] 



(22) 



For a single-mode point contact Ac ~ 1 it is clear 
that the thermodynamic bound is irrelevant. However 
it becomes relevant if one has 20 or more such single- 
mode point-contacts in parallel thermally (for example 
as sketched in Fig. 2b of Ref. (13]). 



XI. NERNST'S UNATTAINABILITY 
PRINCIPLE 

Nernst's unattainability principle, is one of the two 
statements which are together referred to as the third 
law of thermodynamics. It says that no matter how much 
work one does, one cannot take a reservoir initially at fi- 
nite temperature to zero temperature in a finite time (the 



other part of the third law is the claim that entropy van- 
ishes at zero temperature). In the context of a reservoir 
i which is a gas of free electrons in equilibrium, we use 
the textbook result that the heat capacity of a such a 
Fermi gas, Cj = (dQ^/dT) cx T. If the quantum sys- 
tem is sucking heat out over a long time via a DC bias, 
then the fastest rate of heat removal is given by quan- 
tum bound, Jf cx Tf. Assuming this heat flow is weak 
enough that the reservoir remains in equilibrium, we have 



jf = {dQi/dt) = Ci(dTi/dt). Thus 



dt 



-I? 



with C — 1- 



(23) 



Solving this for arbitrary £, one sees that T = is 



reached in a finite time if £ < 1, but not if £ > 1. Eq. ( 23 1 
has £ = 1, which is the critical value of £, then the tem- 
perature decays exponentially with time i; thus it never 
reach absolute zero, but does get exponentially close to 
it. A concrete example with £ = 1 is the point-contact 
considered in Ref. [13] at large bias. For £ > 1 (e.g. the 
pump in Ref. [55] has £ = 3/2 for a fcrmionic bath) the 
decay is much slower, going like 1 > for large t. 

Ref. [30 proposes a quantum system with £ < 1, vio- 
lating the unattainability principle. It is clear that their 
system does not fall into the broad class of systems con- 
sidered here. In particular the bosonic reservoir they cool 
is very different from the ones that I consider. 



XII. CONCLUDING REMARKS 

This work is on the DC response of any quantum device 
where single-electron charging effects are not significant 
(single electron charging energy much less than temper- 
ature or level-broadening due to the coupling to reser- 
voirs). Here "DC" means the zero frequency response 
to time-independent potentials. However, many classical 
heat-engines and refrigerators rely on a pumping cycle, 
so it would be worthwhile to extend this analysis to the 
quantum equivalents of such finite-frequency driven sys- 
tems (quantum pumps, etc). One should also consider 
that quantum transport is inherently noisy, so the energy 
and entropy in the quantum system fluctuates randomly 
on short timescales. These fluctuations give little contri- 
bution to low-frequency responses, because they average 
out on long timescales. In the future we hope to ad- 
dress whether such quantum fluctuations obey the clas- 
sical fluctuation theoreirP4 

Finally, it is worth recalling that all the results pre- 
sented here rely on only two ingredients; (i) Schrodinger's 
equation formulated as a scattering theory, and (ii) sim- 
ple equilibrium statistical mechanics in the reservoirs. 
It is remarkable that the thermodynamic laws for non- 
equilibrium transport (far outside the linear-response 
regime) emerge so naturally from a simple combination 
of these two ingredients. 
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